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В статье представлена моментная схема конечного элемента на примере пространственного 
прямоугольного параллелепипеда. Выведенные вариационные формулы являются статическими 
соотношениями и используются как альтернативный метод расчёта напряжённо-деформи- 
рованного состояния эластомерных конструкций. Они позволяют инженеру-проектировщику 
быстрее И эффективнее осуществлять расчёт исследуемых объектов. 
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Введение 


При использовании традиционных схем метода конечных элементов (МКЭ) 
в форме метода перемещений, построенного на основе вариационного принципа 
Лагранжа для решения задач с особенностями (таких как учёт слабой 
сжимаемости, расчёт пластин на базе трёхмерных конечных элементов и др.) 
возникают существенные трудности [2, 3], для преодоления которых 
используются другие вариационные принципы — Кастильяно (метод сил), 
Хеллингера — Рейсснера, Ху — Вашицу (смешанный метод) и др.[4]. 

МКЭ в форме метода сил не получил значительного развития в силу 
сложности при аппроксимации напряжённого состояния. Большее применение 
нашли смешанные схемы МКЭ. Имея положительные особенности [5], они 
обладают и рядом недостатков, как увеличение порядка разрешающей системы 
уравнений по сравнению с МКЭ в форме метода перемещений, нарушение 
положительной определенности матрицы уравнений. Поэтому для задач с 
указанными особенностями предпочтительнее развитие гибридных схем МКЭ в 
форме метода перемещений на базе вариационного принципа Лагранжа. 

Стандартный МКЭ в форме метода перемещений требует, чтобы поле 
перемещений точек внутри КЭ аппроксимировалось полиномиальными 
функциями, а контакт на границах элементов осуществлялся при соблюдении 
условий неразрывности. Как следует из работ [6; 7] этот вариант МКЭ обладает 
медленной сходимостью в силу полиномиальности функций. Аппроксимирующие 
поля перемещений не включают слагаемое, описывающее жесткие смещения КЭ. 
Этот эффект существенно проявляется при использовании криволинейных КЭ, и 
учет жестких смещений КЭ следует рассматривать не как необходимое условие 
сходимости, а важное средство повышения эффективности МКЭ при расчете тел 
криволинейной формы. 

В процессе эксплуатации стандартной схемы МКЭ в форме метода 
перемещений наряду с проявлением жестких смещений КЭ было замечено и 
другое негативное свойство МЖ, называемое «эффектом ложного сдвига» [1], а 
именно: при изгибе тонких пластин и оболочек на базе трехмерных КЭ 
значительно возрастают погрешности, связанные с проявлением фиктивных 
сдвиговых деформаций. 

Для устранения этих недостатков была разработана моментная схема 
конечного элемента (МСКЭ) [1, 2, 8], позволяющая учесть основные свойства 
жестких смещений для изопараметрических и криволинейных КЭ изотропных 
упругих тел. Суть её заключается в отбрасывании некоторых членов разложения 
деформаций, реагирующих на жесткие смещения и на появляющиеся фиктивные 
сдвиговые деформации. При этом точные уравнения связи деформаций и 
перемещений заменяются приближенными. 

При решении практической задачи механики эластомеров встаёт проблема 
выбора наиболее удачной, по возможности оптимальной, расчётной схемы, 
которая базируется на целом ряде конкретных методов вычислительной 
математики [2]. Но на данном этапе, через недостаточные научные исследования, 
говорить про оптимальность той или другой вычислительной схемы сложно. Это 
обстоятельство часто заставляет строить разные вычислительные алгоритмы, а 
потом сравнивать их преимущества и недостатки. Полученные окончательные и 
некоторые промежуточные результаты расчётов должны исследоваться на 
соответствия их механического смысла задачи. Это так же является необходимой 
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частью расчётов, т.к. ошибки округления и проявление нестойкости некоторых 
расчётных алгоритмов могут значительно изменить результат. Процесс анализа 
результатов есть очень сложным и трудоёмким процессом [6]. 

Выходом из данной ситуации может быть разработка систем 
автоматизированного проектирования (САПР), которое позволило бы инженеру - 
программисту создавать надёжные модели в своей предметной области. Данные 
расчётные комплексы должны состоять из систем прикладных программ и 
содержать в своих библиотеках заранее разработанные модули расчётов 
конкретных объектов [5]. 

Цель этой работы — представление моментной схемы конечных элементов 
для решения задач механики эластомеров и использование для расчётов в 
интелектуальных системах. 


Вывод формулы 


Положим, что область, занимаемая элементом, отображена в кубе с 
единичными рёбрами (рис.1). 


Рис. 1. Шестигранный прямоугольный конечный элемент 


В центр «изопараметрического» КЭ поместим начало местной системы 
координат Оху2, направляя оси вдоль рёбер [1, 4]. 

Функции, задающие геометрию для криволинейного КЭ в базисной системе 
координат, представляется в виде соотношений: 
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где М, = 5а + хх, (1+ уу, (1+ 22,) -функции формы, заданные в базисной 


системе координат; 2",2”, 2“ - координаты узловых точек параллелепипеда. 
Определим якобиан перехода от базисной в местную систему координат: 


дм ОМ, ОМ, 
м м м 
2-м 9 ее 
Эм ОМ, ОМ, 
№ 02 ду | 


(2) 


где 


1- количество узлов; фи 1 абсциссы 1-го узла; ыы } ординаты -го узла; > | 
аппликаты 1-го узла. 

Основной принцип моментной схемы лежит в представлении 
аппроксимирующей функции в ряд Тейлора с последующим отбрасыванием и- ых 
членов ряда. Для пространственного прямоугольного шестигранного конечного 
элемента эта функция, представленная в виде ряда, будет иметь следующий вид 


[2]: 


—_ 000 100 100 010 010 001 001 110 110 101 101 011011 И 
ии=и, чи у чи и чм у Ч у чу чу у 


(4) 


где у,” -коэффициенты разложения, у” - набор степенных координатных 


функций, определяемых по формуле: 


у 


ра и 


раг 


(р=0, 1; а=0, 1;"=0, 1). 
(5) 


Компоненты тензора деформации разложим в ряд Макларена в окрестности 
начала координат: 


р 
— (518), „(58) 
=> е у 


512 


(6) 
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ИЛИ 
8 ‚=е+е и” ‘чел у” ве” 
Ба ное, 
паи рим 
(7) 
поте 
ое нуие 


еб + ео 
83 =@з +ез и 


В разложении компонент о ораняЯ наряду с коэффициентами 
разложения деформаций присутствуют коэффициенты разложения жёстких 
поворотов. Это обстоятельство обуславливает причину замедленной сходимости 
КЭ. Чтобы её устранить отбросим эти члены ряда. После преобразования для 
данного конечного элемента, тензоры деформации будут иметь следующий вид: 


= Уи (и + бу + и 2+ бу); 
1 
(8) 
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11 


8 
г. (ИЯ +, °х) + @, (и +и/:х), 


где и, у, @ - компоненты перемещений каждого узла в пределах КЭ; 


й 
и,- коэффициенты поворота узлов, Ь. „„)“ Коэффициенты, которые связывают 


значения узловых поворотов и степенные функции и (Р9”). Для кубического КЭ, 
который имеет единичные метрики измерения, представленного в естественной 


системе координат, можно определить коэффициенты р . Так, 


Во = Бо = Вы =1, а Ву =Ви = Ви =Ви=0. (9) 


Представим формулу удельной энергии деформации — системы. 
Первоначально преобразуем её к виду: 
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1 1 1 
П = ЛА рьеь +215) +Озёз)+ Е» о Ва аа Ь Рен: 
у 
1 1 1 
+85; ие) +8, Рубь) +8, Рыб) › 


(10) 


где ШО„-коэффициенты матрицы модулей упругости конечного элемента, 
т, п = 1,6, причём О =Ощ. 


Используя (10), получим формулы для вычисления коэффициентов 
матрицы жёсткости [ К |: 


__9?П — 01 „ЗАВ. 

ил ди ди, 2°и]2 ^ дидь, зо ди ди, ) 

_ 6? П 2 ие М 

ЗЛ  дудби,' 3? дубу,” 2 дуди › 
(р 

П П П 


Е.А. ЗАВ д ВА 
"’ диби, “” дб, “” диди, 


Выводы 


Изложенная в статье схема вывода соотношений МКЭ позволяет учесть 
основные свойства жёстких смещений, как для изопараметрических, так и для 
криволинейных конечных элементов, приводящие к таким нежелательным 
последствиям. 

Представленная методика нахождения потенциальной энергии системы на 
основе моментной схемы конечных элементов является универсальной и 
позволяет численно рассчитывать задачи механики эластомерных конструкций. 
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С.И. Гоменюк, В.В. Лаврик 

Представление моментной схемы конечных элементов для расчёта 
напряжённо-деформированного состояния конструкций из эластомерных 
материалов в инструментальных системах 

В данной статье разработаны соотношения модели принятия решения для 
интеллектуальных систем. 

Для реализации данной схемы предлагается совместное использование 
экспертных систем и методов эволюционной адаптации для более эффективного 
решения подобных задач. Показано, что экспертная система позволяют 
разработать структурированную схему, отражающую весь ход процесса принятия 
решений в неопределенных и расплывчатых условиях. Система моделирует и 
интерпретирует действия пользователя по организации своих знаний об объекте и 
делает из них выводы. 

Разработанная схема позволяет совместное использование экспертных 
систем и методов эволюционной адаптации для эффективного решения задач и 
позволяет внедрить разработанную структурированную схему, отражающую весь 
ход процесса принятия решений, в любую интеллектуальную систему. Эта схема 
даст возможность составить конечное множество вопросов, которые помогут 
пользователю эффективно произвести весь процесс программной реализации. 
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